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Abstract 

We analyse the effects that the impermeability constraint induces on the equilibrium 
shapes of a three-dimensional vesicle hosting a rigid inclusion. A given alteration of the inclu- 
sion and/or vesicle parameters leads to shape modifications of different orders of magnitude, 
when applied to permeable or impermeable vesicles. Moreover, the enclosed-volume con- 
straint wrecks the uniqueness of stationary equilibrium shapes, and gives rise to pear-shaped 
or stomatocyte-like vesicles. 
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1 Introduction 



Lipid proteins embedded in biological membranes strongly influence both the geometric and the 
clastic properties of the hosting membranes. The rigid inclusions are able to induce vesicle bud- 
ding m3El> while the interplay between the protein- membrane interaction and the spontaneous 
curvature may yield a loss of regular equilibrium configurations [3]. Moreover, the elasticity of 
the hosting vesicle induces a membrane-mediated interaction that has been widely studied both 
experimentally |S] and theoretically, in the cases of planar |SJ [7], quasi-planar OH E3], and 
quasi-spherical vesicles jllj . 

In this paper we focus attention on the equilibrium shapes of a three-dimensional vesicle 
hosting a rigid inclusion, and in particular on the effects of the permeability properties of the 
vesicle. In fact, a given slight perturbation applied to a quasi-spherical vesicle may induce quite 
different changes on the resulting equilibrium shape, depending on whether the enclosed volume 
of the vesicle is constrained or not. More precisely, an O(e) relative perturbation of the external 
parameters induces an equivalent O(e) modification in the shape function of a permeable vesicle, 
but a quite stronger 0(y/e) relative perturbation if the enclosed volume is kept fixed. Furthermore, 
the volume constraint yields a multiplicity of stationary equilibrium shapes, and it leads to abandon 
the spherical shape towards either pear-shaped or stomatocyte-like vesicles. 

The volume enclosed by the vesicle may be constrained or not depending on both the chemical 
properties of the aqueous solution which surrounds the vesicle, and the time scales in which the 
(meta)equilibrium shapes are observed ^JEH- When the water is essentially free of molecules 
that cannot permeate the bilayer membrane, no volume constraint stands. On the contrary, when 
some of the molecules in the solution are unable to permeate the bilayer, the resulting osmotic 
pressure gives rise to an enclosed-volume constraint. However, even in this latter case, on long 
time scales water molecules succeed in permeating the membrane. Eventually, the vesicle reaches 
its true equilibrium shape, which minimizes the free energy with respect to the enclosed volume. 

Throughout our development, we will model proteins as rigid inclusions, and we will assume 
that the protein- membrane interaction simply fixes the contact angle, i.e. the angle the vesicle 
normal determines with the inclusion plane |14l l%llTK ] . However, only slight modifications need to 
be applied to our results to take into account interactions which determine the contact curvature 
instead of the contact angle ^JEl^j- More drastic changes, though substantially the same 
mathematical setting, requires the weak anchoring case ^3 ^3 H i m which the inclusion- vesicle 
interaction provides an additional term in the free-energy functional, instead of a fixed boundary 
condition. 

We describe vesicle elasticity through Helfrich's (spontaneous curvature) model |2()l I21| . The 
free-energy functional is then 

:=k [ (H-a ) 2 da , (1.1) 

where £ is a closed surface describing the vesicle shape, H denotes the mean curvature along £, 
k is the bending energy, and ag the spontaneous curvature. In the minimizing procedure at fixed 
area (and possibly fixed volume) we replace (|1.1|) with the effective free-energy functional 

.Fes [£] := k / (H - a ) 2 da + X (Area {T.) - A) + L (Vol (£) - V) | . (1.2) 

The Lagrange multipliers A, fx have the physical meaning of surface tension and pressure difference, 
and the brackets are there to remind that the volume constraint will not be always applied. In 
(|1.2[l the area constraint has been inserted as a global, instead of a local constraint. We recall 
that, in the absence of external forces, both choices are equivalent |22j . 

The Euler-Lagrange equation associated to the functional (|1.2|) is the shape equation [221 123] : 

K [A S H + 2H(H 2 -K) + 2a a K - 2a^H] -2XH- L| = , (1.3) 

where A s , the tangential divergence of the tangential gradient, is the Laplace-Beltrami operator 
on E, and K denotes the Gaussian curvature along S. 
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Figure 1: Cross-section of a three-dimensional vesicle embedding a conical inclusion of base radius 
a. The tangent to the generating curve at the contact points determines a fixed angle ip with the 
inclusion plane. 



The plan of the paper is as follows. In the next section, we introduce the surface parameteri- 
zation and we derive the conditions satisfied by spherical and quasi-spherical shapes. In Section 3 
we analyse the vesicle shapes obtained in the presence of the area constraint alone (i.e., the long- 
time equilibrium vesicle shapes). Section 4 is devoted to the peculiar role played by the volume 
constraint. In Section 5 we review and discuss our results. 



2 The model 

Let us consider a vesicle which embeds an inclusion, that we model as a symmetric conical frustum 
of negligible height, base radius a, and apex angle tp. The inclusion-vesicle interaction fixes the 
angle between the vesicle normal and the inclusion plane at the contact points to be equal to 
(f ~~ v)' Figure ^ illustrates the geometric setup of the model. For a more detailed description 
of the inclusion- vesicle interactions and their modeling we refer the reader to the paper by Biscari 
and Rosso @|. 

We restrict our analysis to axisymmetric vesicle shapes, and parameterize them in spherical 
coordinates centered in a fixed point O, which lies in the inclusion symmetry axis z (we assume 
that this is the symmetry axis of the vesicle too, since the inclusion does not upset the cylindrical 
symmetry) : 

P("&, if) — O = r(i?) sin d cos f e x + r(i?) sin d sin ipe y + r(i?) cos $ e z . 

The unit vectors {e^, e yi e z } form an orthogonal basis, with e z parallel to the inclusion axis; if) 
are the polar and azimuthal angles, respectively. The area element, the Laplace-Beltrami operator, 
the mean curvature, and the Gaussian curvature along £ are given by: 



da = J~g did dip , A s = — 

V9 



r sin $ \ / \J r 2 + 
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2r 3 + 3rr' 2 - r' 3 cot ■d - r 2 (r' cot + r" ) 
H = ttt; , and 



K = 



2r (r 2 +r' 2 ) 3/2 
(r sin d-r' cos ■&) (r 2 + 2r' 2 -rr") 
r sin d (r 2 + r' 2 ) 2 



where a prime denotes differentiation with respect to the polar angle, and g := (r 2 +r' 2 ^)r 2 sin 2 •& 
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The spherical parameterization transforms the effective free energy <|1.2|) as follows: 



F cS [ r ] = 2ir / (re (H[r] - cr ) + A) r y/r 2 + r' 2 sintfdi? + — £t / r 3 (i?) sin eft? 

The integration limit #f depends on position of the origin O. We choose to fix O at a distance 
Az = (a cot ?/0 above the inclusion, which yields $f = (n—ip). This choice simplifies the attachment 
condition on the inclusion and the constraint on the direction of the contact normal, which become: 

r (i9 f ) = and r'O&t) = . (2.1) 

smtp 

The vesicle is free on its top (i? = 0). The boundary conditions therein follow from regularity 
requirements on the vesicle shape: 

lira r'(i?) = and lim r"'(i?) = . (2.2) 
2.1 Internal actions 

Let us consider a subsurface £' C E, and let n and t respectively denote the normal to E at a 
point P 6 $£', and the tangent to the curve <9£'. Furthermore, let k := nAt denote the direction 
in the tangent plane at P pointing outwards with respect to £'. The internal actions at P consist 
in a distributed force f and a distributed torque m, whose densities per unit length of <9£' are 



k(H-ctq) 2 + A 



dH 

k — — n and m = n (H — erg) n 
ok v ' 



The above equations generalize to three-dimensional vesicles the internal actions derived in |7] in 
the two-dimensional case. They show that the surface tension A may become negative without 
giving rise to the collapse of the vesicle, provided that f k = [re (H— CTo) 2 + A] remains non-negative 
all along the vesicle. Furthermore, they provide an alternative way of deriving the free-boundary 
conditions 12. 2|) . In fact, these conditions arc equivalent to the vanishing of the internal force and 
torque acting on an infinitesimal cap which surrounds the vesicle top. 

2.2 Spherical shapes 

The equilibrium shape of a vesicle is a sphere of radius r , centered in O, whenever the base radius 
a , the apex angle tp , and the vesicle area A a satisfy 

clq = t sinV'o and A Q — Airr^ cos 2 — . (2-3) 
If we eliminate r Q from equations Q2.3[l i and l|2.3|) 9. we obtain: 

2 a ■ 2 ro 
ira — A sm — . 

If, in addition, the vesicle is impermeable, the enclosed volume V a must match 

4 1 2 1 9 , TOo COS% f . 

Vo = -Krt cos — + -ira r cosip = — , g (2 - cos^oj. 
3 2 3 6 S m J ^ 

In particular, area, enclosed volume, and apex angle must obey: 

wo := 36tt-§ = (2 - cos ifj ) 2 cos 2 ^ . (2.4) 
A% 2 

We remark that for any value in v 6 [1, 2] (the end cases corresponding to the cases of a sphere 
and a half-sphere), there is exactly one value of ip € [0, 5] that satisfies H2.4|) . 
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2.3 Quasi-spherical shapes 

We now assume that some of the control parameters a, ip, A (and possibly V) are slightly perturbed 
with respect to their values satisfying (|2.3|) and l|2.4l) . In this case, we look for solutions of the 
shape equation (|1.3|1 that represent a perturbation of a sphere: 

r(0) =r o (l + e$i(0)+o(e)) . (2.5) 

Consistently, we also expand the Lagrange multipliers by perturbing their "spherical" values: 

A = 4 (Ao + eAi +o(e)) and /i = ^ (770 + e 771 + o(e)) . (2.6) 

In H2.6fl . the normalizing factors n and r Q have been inserted in order to proceed with the dimen- 
sionless quantities Aj and r\i. Furthermore, in order to make the whole shape equation dimension- 
less, we define the reduced spontaneous curvature as 

:= <to r ■ (2.7) 
If we insert 1)2. 5p. (|2.(i[) and (|2.7|l in \1.'S\ . we derive at O(l) the condition 

A + »»=«)(l-ft) ) (2.8) 

linking the spherical values of the Lagrange multipliers. 

When we push the expansion to O(e), we obtain a fourth-order linear differential equation for 
Qi. If we further introduce the variable s := cost?, and perform the substitution gi(t9) = n7i(cosi9), 
the so-obtained differential equation reads as: 

(1 - s 2 ) 2 n7 (4) - 8s (1 - s 2 ) zu {3) + [12s 2 - 4 + 51 (l - s 2 )] w {2) - 2g 1 szu < - 1) + g a zu = -4(A X + 771), 

(2.9) 

where 

30 = 251 := 4(<to(2 - cro) - A ) , 

and the superscripts denote differentiation with respect to s. Equation (|2.9() is an inhomogeneous 
fourth-order Legendre differential equation. Its general solution can be expressed in terms of 
Legendre functions of the first and second kind as: 

= _ 4(A 1 +r ?1 ) + ^ + ^ + Ca p ^ + ^ Q ^ 

50 

where the orders of the Legendre functions are given by 



1 1 



v± = - g + 2 V 5 + 2 -9i ± + 2) 2 - 4 5 o • (2.10) 

1 1 



If we replace go — 2gi in Ij2.1t)|l . we obtain u± = — — + — y5 + 2gi ± 2 |<?i — 2| , so that 



{^+^-} = {i^(v / T+4ffT-i)} 



If we further introduce the notation 

v := ~ (y/T+isn- l) = i f v/1 + 8 ? o(2 - to) - 8A - l) , (2.11) 



and we finally consider that -Pi(s) = s , we arrive at the general solution of the linearized shape 
equation: 

ei (A) = _ 4 (Ai +m) coB^ + cijQ^coa^ + CaP^rcos^ + GiQvfcosi?) . (2.12) 
.9o 
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Some further investigation will turn out to be necessary in the particular cases v = (i.e. g\ = 0) 
and v = 1 (i.e. g\ — 2). We postpone the analysis of these cases to the sections below. 

In the following, we will analyse and compare the perturbed equilibrium vesicle shapes of 
permeable and impermeable vesicles. The derivation below works when any or even all of the 
parameters a, ip, A, or V are varied with respect to their spherical values. However, and only 
in order to shorten our presentation, we will henceforth restrict our development to the case 
in which only the area, and possibly the enclosed volume, are varied with respect to the values 
satisfying l|2.3(l and (|2.4(1 . while the inclusion parameters are kept unchanged. These are the easiest 
perturbations to be implemented experimentally: for example, an area variation in a biological 
membrane may be simply induced by adding extra lipid molecules to the membrane bilayer. 



3 Permeable vesicles 

When the vesicle is inextensible but permeable, no volume constraint stands. For all practical 
purposes, this is equivalent to assume that the Lagrange multiplier // vanishes identically. When 
this is the case, condition (|2.8|) reads as 

A = ?o(l - Co) , 

which implies go — 2g\ — 4c , and 

v = \{Vl + %^-l) • (3.1) 
The linear perturbation of the spherical shape becomes: 

Qiiti) = -— +Ci cost? + C* 2 Qi(costf) + C* 3 P„(costf) +C 4 Q v {cosd). (3.2) 

All Legendre functions of the second kind Q v {s) are singular when s — > 1~. Thus, the free- 
boundary conditions l|2.2|) require C2 = C4 = 0. The remaining parameters Ai, C\, and C3 can 
be determined with the aid of the contact conditions l|2.1|l and the area constraint: 

Ql ($ { ) = 0, ^(^=0, and 4nr% qi (t?) sin dti = AA , (3.3) 

Jo 

where A A is the area excess with respect to the spherical value A . With the aid of l|A1.2|) - (|A1.4|l 
conditions l|3.3fl yield 

Ai = Co C3 esc 2 ip [V cos^o P v -x{— cosijj ) + (sin 2 %p a + vcos 2 ip Q ) P v {— cosV'o)] , 

Ci = — vCz esc 2 tp \P v -i{— costAo) + cos ip Pv(— costAo)] , and 



C a 



4(P I/ _i(- cosi/>o) - P v {— COSl/'o)) 



(v + 1)(1 — COS tpo) 

— {v — \){v cos V^o + 2)P V {— cos?/> ) — [y A + v + 2)P v -i{— cosipo) 



1 + COS tfjo 

AA 



A c 



We have already announced that the differential equation (|2.9J) admits l|3.2(l as its general 
solution only when v ^ {0, 1}, that is when the spontaneous curvature is neither null nor equal to 
the inverse of the unperturbed radius r . We will now solve l|2.9|l in these cases. 

When <;o = 0, the general solution of the homogeneous differential equation associated to (I2.9|l 
is still as in (|3.2() . with v — 0. However, and since Po(s) = 1, the particular solution of the 
equation is not a constant. By using the method of variation of parameters, and requiring also 
the free-boundary conditions H2.2(l . we find: 

g (w=0) = Ci cos ^ + C 3 + Ai (3 + 2 lg(l + cos )) . 
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In addition, conditions l|3.3|) yield: 



2 "00 



Ai = -d sin 2 ^ , C 3 = Cx 



2 V'o \ \ . 2 V'o 



21g(2sin 2 ^ I I sin , ; 



2 i>° , o 4. 2 V'o , . 2 \ ^ AA 



and 

2(1 + sin 2 — + 2tan 2 — lgsin 2 — ) d = ^ 
V 2 2 5 2 / A 

When = 1, the particular solution of the differential equation (|2.9(l is again a constant, but 
the solution of the homogeneous equation is not of the form (|3.2I) . If we solve l|2.9[) explicitly and 
use the free-boundary conditions l|2.2[) . we arrive at: 

^ 0=1 > (tf) = Ci cos d + C* 3 ( cos i? lg(l + cos 0) - 1) - Ai . 
Finally, conditions 1)3.3(1 now require: 

a 1 = -C 3 flg(l-cos^)-T^V > ) - Ai=-^fl + T^r > ) ■ and 



1 — cos tp J ' \ 1 - COS Ipo 

2 ,2 Ipo n AA 

cos T cot Y C 3 = — . 

Figure shows how the perturbed vesicle shape depends on the spontaneous curvature for a 
prescribed area increase (5% with respect to the spherical value), while Figure shows the volume 
variation induced by AA. 

• Since the inclusion is placed at § = i?f , Figure [3 shows that in vesicles characterized by 
greater spontaneous curvatures the shape modifications gather away from the inclusion. 

• An area increase induces a shape perturbation g± that does not change sign all along the 
vesicle (we will find below that this is not the case when also the volume is constrained). 

• Figure |3| shows that the volume increase induced by AA increases monotonically with AA. 
However, the remarkable increase in AV^ that shows up when Q) ~ 3 is to be linked with 
the spontaneous-curvature driven budding transition |24j that is close to occur. A detailed 
analysis of the inclusion's influence on this transition can be performed only by studying the 
nonlinear shape equation 1)1. 3J1 . and will be reported elsewhere |25| . 



3.1 Small inclusions 

By using the asymptotic expansion (jAl . it is possible to show that in the small inclusion 
limit a <C \/A2, which implies ipo "C 1 by virtue of 1(2. 3|) , the perturbed vesicle shape becomes 
independent of 

Figure [5] shows that when the spontaneous curvature is small, the asymptotic expression ((3.4(1 
is more rapidly approached. In fact, if we compute the volume variation associated with the 
perturbed shape 1)3.4(1 we obtain: 

AF 3 AA / (2l , fAA\ 2 \ 

— 2x +0 (^ lg Mx) J • 

Figure |31 confirms that the relative volume variation is closer to the small limit prediction ^AA/A 
when is small. 



G 



2ZL &{ 



Figure 2: Perturbed shapes of a vesicle hosting an inclusion, when the area is slightly greater than 
the area corresponding to a spherical equilibrium solution. The plots correspond to tp — O.I71", 
AA = 0.05 A , and Q) = 0, 1, 2, 3 (the arrow points towards increasing values of Q)). 



.080 



AV/V 



.079 




Figure 3: Volume variation induced by a given area variation in a permeable vesicle hosting an 
inclusion, as a function of the reduced spontaneous curvature Co- As in Figure [21 1^0 = 0.17T and 
AA = 0.05 A a . 
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4 Impermeable vesicles 



We now focus on the solutions of the differential equation (|2.9[1 that satisfy both area and volume 
constraints, when these geometrical quantities are close to satisfy the spherical condition 12. 4|) . 
Leaving aside for the moment the cases v = and v = 1 (that turn out to be meaningless in the 
impermeable case), the solution of l|2.9[) is of the form 12.12|l with the following parameters to be 
determined: C\, C2, C3, C4, Ao (i.e. v : by virtue of l|2.11|l ) . and the combination (Ai + 771). 



4.1 Singular perturbations 

The boundary conditions 1)2. lfl and 12.211 may determine the four parameters C1-C4. A problem 
arises when we try to determine v and (Ai + by using the area and volume constraints. In 
fact, the 0(e) of the above constraints reads as: 

i5f A A r#t \y 

gi (d)smftdft = - — / gi{ft)smftdft = -—^ . (4.1) 
It is clearly impossible to satisfy both constraints if 

AA^—AV. 

To 

From the analytical point of view, the degeneracy of the constraints (|4.1|l stems from the fact that 
the first variations of area and volume of a surface are linearly dependent when computed in a 
spherical shape. In fact, if we perturb any surface with constant mean curvature H (the so-called 
Delaunay surfaces |26l 127) ). the area and volume of the resulting surface satisfy [231 

AA = 2HAV. 

The area and volume constraints become linearly independent only when the second variations 
come into play. Thus, if we are willing to perturb the assigned area and volume to an arbitrary 
0(e), we have to perturb the shape function r(ft) to 0[y/e), as we shall show next. 

We begin by replacing (|2.5II by 

r(ft) = r (l + ^0i(tf) + e£i(tf)+o(e)) . (4.2) 

The singular perturbation gi is of the form l)2.12[l . with Ai and r\\ replaced by their half-order 
counterparts A 1 and rji . Correspondingly, the area and enclosed volume of the resulting vesicle 
shape are given by: 

A + eAA = 2n r \J r" 1 + r' 2 sin d dd 
Jo 

= A a + 4nrly/e J qi sin ft dd + -Kr 2 a e J (2q\ + gf + 4g^j sinftdti + o(e); (4.3) 

V +eAV=— r 3 (tf)sintfdtf = 

= V + 2Trrl^J gismi9d$ + 7rr 3 ej^ (2g | + 2g^ sin t? d<& + o(e) . (4.4) 
Both H4.3f) and (|4.4|l can now be satisfied, provided the functions gi and Q\ are such that: 

gi(<&) sintfcii? = , (4.5) 



g'!W-2gl(d) sm$d$= To ^ , and (4.6) 



AA - 2AV 



#t AT/ 1-0 



AV 1 

gd-d) smft dft = — - / gl(ft)svaftdft . (4.7) 
2irr-l 



S 



Equations (|4.5|) and 1(4.6(1 fix g 1 , as we will show below; then, equation (|4.7J) allows to determine 
also the now next-order correction g\ . 

Before entering in the detailed analysis of equations (|4.5|l and 1(4. 6|1 . the above result deserves 
some remarks. 

• Equation 14.2(1 underlines the most striking effect of the impermeability constraint on the 
vesicle: an O(e) perturbation of the assigned geometrical values induces an (9(v / e) pertur- 
bation in the vesicle shape. 

• Equation l|4.5|l is an eigenvalue equation. We have to look for non-trivial (i.e., non- vanishing) 
perturbations gi that satisfy it. It will prove to admit a countable infinity of independent 
solutions. 

• Equation 1(4.61) is quadratic in g 1 . We will thus find two, rather than one, possible perturbed 
shapes for any non-trivial solution of (|4.5|) . An energy argument will be needed to identify 
the preferred perturbation among the double-infinity of possible choices at our disposal. Both 
perturbed shapes arising from 1(4.6(1 will deserve notice, since they display two qualitatively 
different vesicle reactions to the perturbation. 

• The quadratic expression in the left-hand side of (|4.6|l forces the right-hand side combination 
(r AA ~ 2AV) to assume only non-negative values. This property is not peculiar of vesicle 
theory: it reflects a classical isoperimetric inequality. In fact, for any closed surface, the 
ratio A 3 /V 2 is bounded from below by the value 36tt, which is attained only by a sphere 
(see, e.g., [25]. p. 8). Thus, for example, it does not exist a closed surface with the same 
enclosed volume of a sphere and smaller area. 

4.2 Multiplicity of stationary perturbed shapes 

The free boundary conditions ((2.21) require that the coefficients of the singular Legendre functions 
of the second kind in (|2.12l) must be null: C-z = C± = 0. Afterwards, the contact conditions 1(2. 1() . 
which are linear in the shape function, supply two relations that connect C\, C3, and (Ai +771) 
in H2.12fl . As a result, the leading perturbation gi in 14.2(1 can be given the form: 

pi (0) = C 3 [a(u,ip ) +/3(^,Vo)costf + P„(costf)] , with (4.8) 
a(v,ip ) := — esc 2 tp (i/cos^ P v -i{— cos ip ) + (sin 2 tp + vcos 2 ipo)Pu{— cos V'o)) and 
0(v, ip ) := -i/csc 2 lp (P !/ _i(- COS^q) + cosV'o P v {— cos %jj )) . 

In order to determine completely the function gi, we have to use 1(4. 5|l and ((4.6(1 to derive v and 
C3. Any non-trivial solution of ((4.5(1 possesses C3 7^ 0. Thus, we can drop C3 out from it, to obtain 
an eigenvalue equation in depending only on ip a . Figure ^illustrates the numerical solutions of 
((4.5(1 . These solutions exhibit the following properties: 

• The spontaneous curvature does not enter in ((4.5(1 . Thus, the stationary shape modifications 
of an impermeable vesicle do not depend on <^o, as they did in the permeable case (see ((3.1(1 
and Figure [21 ■ 

• For any ip G [0, there is a countable infinity of values k%(^o) satisfying ((4.5(1 . 

• For any ip £ [0, -|], the solutions ^fc(V'o) are symmetric with respect to v = — |. However, 
since the symmetric solutions are identical (see (|Al.l(l 'l. we can restrict our attention to 
solutions with v > — ^. 

• When f <l (small protein limit), the use of ((A1.5(l allows to prove that tends to an 
integer value for any k: 

VkU>o) = (k+l) + £^%2 + o(V2) for any fceN. (4.9) 
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The shape function thus approaches a linear combination of Legendre functions of integer 
order, that is, Legendre polynomials. Furthermore, both l|4.9|l and Figure 0] show that only 
Legendre polynomials of order equal to or greater than 2 come into play, while only low- 
order Legendre polynomials (Pq and Pi) entered the small-protein limit of permeable vesicles 
(see 113.411 s ). This yields more drastic shape modifications in the incompressible case, since 
Legendre polynomials are more and more oscillating as their order increases. The expansion 
of the shape function (in the absence of inclusions) in terms of Legendre polynomials was 
first used in |24| . 

• For any k £ N, the functions fki^o) increase monotonically with tp , and do not intersect. 

• It is possible to prove by direct inspection that v — and v = 1 do not solve 14.5fl for any 
value of ip . 




Figure 4: Order v of the Legendre functions entering in the perturbation of an impermeable vesicle 
shape, as a function of the inclusion apex angle ip . The graphs display the smallest five numerical 
solutions {z/fc, k = 1, . . . , 5} of the equation l|4.5fl . 

Once we have identified the z;- values that satisfy l|4.5[) . we can insert (|4.8|l in (|4.6|l to determine 
C3, the only remaining free parameter in the singular perturbation gi. We stress again that, 
being gi linear in C3, this latter parameter enters quadratically in (|4.6(l . This fact, on the one 
hand fixes a sign for the geometrical quantity (r AA — 2AV), and on the other hand implies 
that, for any k S N with v = Vk(jjj ), there are exactly two values of C3 (one the opposite of 
the other) that satisfy 14.6(1 . Thus, for any k E N, there are two possible perturbed shapes: 

r k± {d)=r (l±Vif>i, fe W + 0(e)). 

Only an energy estimate can help us in determining, for any tj) S [0, both the value of v k 
and the sign of C3 that minimize the elastic energy. This will be the aim of the remaining part of 
this section. 

4.3 Energy estimates 
4.3.1 Ground state energy 

In order to identify which is energetically preferred among the stationary perturbed shapes deter- 
mined above, we will now compute their elastic energy. When we insert (I4.2|l in the free-energy 
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functional and we make use of (|4.5|) - (|4.7|) . we obtain: 



[H k ± - (Jo)" r k ± \/rl ± + r'£ ± sintfcW = 



k(to - 1) ^2- + k?o(«> - 1) — - «<To - 3 + F£ + o(r AA, AV) 



AA 



r a AA - 2AV 



where T k ± denotes the free energy of the fc-th solution of equation 1)4. 5[l . with the positive or 
negative sign for C3, and 



•Fl ■= -z K 
2 



sin i9 di9 



(2) 

is the leading order that may determine k and the sign of C3. However, Tj, depends quadratically 
on gi k , so that the sign of C3 cancels out from it. Thus, the second order expansion turns out 
to be able to identify only the preferred value of v k ■ A further term in the free-energy expansion 
will be needed to complete the determination of the free-energy absolute minimizer. 

Figure shows how depends on k and ip - the free energy increases when either of these 
increase. In particular, Figure proves that the stable perturbed shape for a quasi-spherical 
impermeable vesicle corresponds to the solution with the smallest possible order of the Legendre 
functions. This is not surprising from the physical point of view, since Legendre functions wrinkle 
when their order increases, and these oscillations increase the elastic energy. 



(4.10) 
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Figure 5: Results of the numerical computation of the integral in (|4.1()|l . when gi k is given in 
(|4.8|) . and v is the fc-th solution of equation (|4.5|) . The graphs display the results for k = 1, . . . , 5. 



4.3.2 Pear-shaped or stomatocytes? 

We still have to choose the preferred sign for the parameter C3. Figure [S] shows the deep, quali- 
tative, differences between permeable (a), and impermeable ((b) and (c)) stationary shapes that 
arise from the same parameter changes: AA — jqA q for all the shapes; AV is left free in (a), 
while it is kept null in (b) and (c). Permeable vesicles modify their enclosed volume in order to 
keep an almost-spherical shape. On the contrary, impermeable vesicles move towards pear-shaped 
or stomatocyte-like equilibrium shapes |29j , depending on the sign of C3 . 

In order to compare the free energies of pear-shaped and stomatocyte-like vesicles, we need to 
push further our free energy expansion. The next order depends on the third power of g± and its 
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(a) (b) (c) 

Figure 6: Perturbed shapes for a vesicle embedding an inclusion of negligible size [ip <C 1) when 
a 10% increase in the vesicle area is imposed (with respect to the value leading to a spherical 
shape). Picture (a) shows the equilibrium shape of a permeable vesicle, which is allowed to adapt 
its enclosed volume. Pictures (b) and (c) refer to an impermeable vesicle when the positive or 
negative sign for the parameter C3 is chosen when solving equation (|4.6II . The inclusion (not 
visible in the pictures) sits in the bottom end of the vesicle. 



derivatives: it is 0(r AA, AV) 3 ^ 2 . By making use of the constraint requirements (|4.5() - (|4.7|l . it is 
possible to obtain: 

•W) = «(« - I) 2 £ + (* - 1) ^ - «b r ° AA ~ 2AV + ^ + + o(r Q AA, AV) 3/2 , 

' o ' o ' o 

where all the terms up to J 7 ^ do not depend on the pear shaped vs. stomatocyte choice, and 

^^( ^^^ (/W+^W) • (4.11) 

In ^4. lip , the plus sign corresponds to the pear-shaped vesicle, the minus sign describes a stoma- 
tocyte, and ?o denotes as usual the reduced spontaneous curvature. Plots of the functions / and g 
(whose explicit expressions can be found in Appendix A2, equations HA2.2MA2.3jl 1 are snown m 
Figure [7| Both / and g are negative for all values of rp - Thus, for any positive value of the spon- 
taneous curvature, the pear-shaped vesicle is the absolute minimizer of the free-energy, whereas 
the stomatocyte represents only a relative minimum. A transition between the two stationary 
shapes can be observed only when negative spontaneous curvatures are induced; more precisely, 
the stomatocyte-like phase is preferred if 

T0<-4t4= : ^cr(V'o). (4.12) 

9 Wo) 



Figure also shows how the critical value of the reduced spontaneous curvature depends on the 

3 

5 



apex angle. In particular, — | < <ro,cr(V ! o) < — f for all values of ip. 



5 Concluding remarks 

In this paper we have analysed how the impermeability constraint may induce singular stationary 
vesicle shapes, and how the embedding of an inclusion may vary the vesicle topology, promoting 
pear-shaped geometries that anticipate critical phenomena such as budding or vesiculation |3(JU31| . 

Our analysis has been based on the linearization of the shape equation close to a spherical 
shape. This approximation does not allow to approach the aforementioned transitions, but in 
turn yields analytical results which prove that the vesicle reaction to a variation of the external 
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Figure 7: Plots of the functions / (dashed) , g (dotted line) , introduced in (|4.11|l , and the critical 
spontaneous curvature Q),cr (continuous line), defined in (|4.12[l . all as functions of the apex angle. 

parameters may not be analytical. In Section 4 we have shown that an O(e) relative variation in 
the vesicle area may induce either an O(e) or an 0(y / e) relative variation in the shape function, 
depending on the permeability properties of the vesicle and the aqueous solution that surrounds 
it. 

The present analytical study is being currently completed by a numerical study |25| . which 
detects how the presence of an embedded inclusion modifies the phase diagrams describing the 
vesicle topology, and under which conditions it anticipates budding phenomena. 
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Appendix Al: Properties of Legendre functions 

The Legendre functions of the first and second kind (respectively denoted as P v (s) and Q v (s)) are 
the solutions of the linear differential equation 

(1 - s 2 ) y"{s) - 2s y'(s) + v{v + 1) y(s) = . 

The Legendre functions of the first kind are regular when s — ► 1~ (with P v (l) = 1), while the 
Legendre functions of the second kind are singular close to both s = ±1. The properties we use 
in this paper are the following (see §8.2, and §8.5) 
For any i/fE and s e [—1,1], 



2^ + 1 



■P«-l{8) W^-l) 



u+1 v ' v + 1 
(l-s 2 )P^(s) =vP v - l {s)-vsP v {s) . 
Equations l|A1.2() and l|A1.3() imply 

-P^-^cos??) — cos $ P v (cos ■&) 



(Al.l) 
(A1.2) 

(A1.3) 



P^cosi?) sintfefr? = - / P v {s)ds = 



v + l 



(i/^-l). (A1.4) 
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For any v E R + \ N, the Legendre functions of the first kind admit the following asymptotic 
expansion: 

T[V + 1) 

where T, "J*, and 7 respectively denote the Euler gamma function, the digamma function, and 
Euler's constant. 



Appendix A2: Third order expansion of the free energy 

The derivation of the third-order term in the free-energy expansion requires a third-order expansion 
of the shape function. Thus, (|4.2|l has to be replaced by: 

However, the higher-order terms Q\ and g$ turn out to enter in the free-energy expansion only 
through combinations that can be related to integrals of gi and its derivatives by making use of 
the area and volume constraints, as it already happens in the second-order expansion (see l|4.7ll ). 
More precisely, the third-order expansion of the constraints yields 



23(tf)sintfdtf + 2 / gi {tf) g± (■&) smtfdti 



1 

10 Jo 3 Jo 

By using (|A2.1|I it is long but straightforward to prove that 



(i?)sim?di? 



(A2.1) 



A 



AA 



^(ps,Bt) = - 1) — + k<; (^o ~ 1) — - «?o 



,AA - 2AV 



T {2) + T {3) + o(r AA, AV) 



.3/2 



with 
where 



•S lrt )=^o) + ftG(^), 



G{ip ) = 2ir 



v k (l + v k ) 3 

n Sk 



COS 2$ , 9 COt ,q 

e i p i + — ^— 



l-cos2i? 5 2 ' 6 



2 5 2 



The functions /, g introduced in l|4.11|l are related to F, G through 

'r AA- 2AV\~^ 



/(^o) = Fty ) sgn(C 3 ) 



/( f v. = (V| J v.^n(r :i ) 1 2AF ) 



sin 1? dd 



(A2.2) 
(A2.3) 
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